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The  complex  ray-tracing  theory  developed  earlier  [J.  Opt.  Soc.  Am.  A  5,  200  (1988)]  is  applied  to  a  study  of  the 
structure  of  the  electromagnetic  focal  field  spectrum  of  paraboloidal  mirrors.  This  problem  was  treated  recently  by 
Barakat  (Appl.  Opt.  26,  3790  (1987)],  who  used  the  Gaussian  vectorial  diffraction  method.  In  this  paper  a  more 
general  solution  to  the  problem  is  presented  that  uses  the  Stratton-Chu-Silver  integral.  The  reflecting  kernel  then 
explicitly  incorporates  the  surface  function,  the  surface  physical  parameters  and  their  dependence  on  the  state  of 
polarization  of  the  incident  radiation,  and  an  unrestricted  aberration  function.  Numerical  results  are  presented  for 
mirrors  made  of  four  different  metals,  and  their  use  in  optimizing  the  mirror  design  is  demonstrated. 


INTRODUCTION 

In  part  I  of  this  series  a  complex  ray-tracing  theory  was 
developed  for  the  solution  of  diffraction  problems  of  com¬ 
plex  conducting  surfaces. 1  In  part  II  the  theory  was  applied 
to  the  ellipsoidal  scatterers.2  Recently  Barakat3  studied  the 
structure  of  the  focal  spectrum  of  a  paraboloidal  mirror, 
assuming  that  the  mirror  is  absolutely  reflecting,  the  inci¬ 
dent  rays  are  symmetrical  and  meridional,  and  the  system  is 
aberration  free.  In  what  follows  it  is  shown  that  these  as¬ 
sumptions  can  easily  be  eliminated  without  adding  much 
complexity  to  the  method  of  solution. 

First,  if  the  mirror  is  metallic,  the  real  refractive  index 
varies  with  the  angle  of  incidence.  For  example,  at  radiation 
wavelengths  X  =  0.1  am  and  X  =  1.315  am  the  real  refractive 
index  of  aluminum  is,  respectively,  14.64  and  1.00  at  normal 
incidence  and  14.66  and  1.06  at  grazing  incidence.2  Thus 
the  reflectivity  and  absorptivity  vary  with  the  angle  of  ray 
incidence  upon  a  metallic  surface,  and  both  vary  with  the 
state  of  wave  polarization  for  both  metallic  and  dielectric 
surfaces.2-4  Therefore  the  conservation  of  energy  in  a  ray 
tube  [Ref.  5,  Eq.  (2.11)]  and  that  carried  out  by  Barakat  [Ref. 
3,  Eq.  (14)]  should  yield  the  polarization  and  angular  depen¬ 
dence  of  the  surface  reflectivity  in  the  aperture  function, 
q (8).  Second,  asymmetrical  reflective  focusing  is  of  more 
interest  in  practice.  Inertial  confinement  fusion  devices,  in 
which  the  laser  beams  strike  the  focusing  mirrors  obliquely, 
and  astrophysical  observatories,  for  which  unobstructed  im¬ 
ages  are  desired,  are  just  two  practical  examples.  Because  in 
Barakat's  approach  the  Kirchhoff  integral  is  evaluated  on  a 
Gaussian  sphere  centered  at  the  Gaussian  focus,  asymmetri¬ 
cal  focusing  can  hardly  be  analyzed  by  his  approach.  Third, 
the  intended  inclusion  of  some  aberration  orders,  third-or¬ 
der  spherical  and  third-order  coma,  can  easily  be  replaced  by 
the  whole  aberration  function  by  using  the  Stratton-Chu- 
Silver  integral  in  place  of  the  Kirchhoff  integral. 

In  this  paper  a  general  treatment  is  presented  for  the 
focusing  of  a  plane  wave  and  a  Gaussian  laser  beam  by  an 
aluminum-coated  paraboloidal  mirror.  Three  other  metal 

0740-3232/88/091 444-06$02.00 


coatings  are  studied  for  comparison  with  aluminum.  The 
surface  geometrical  and  physical  parameters,  the  state  of 
radiation  polarization,  and  the  whole  wave-front-aberration 
function  are  explicitly  incorporated  into  the  reflecting  ker 


MATHEMATICAL  FORMULATION 

The  equations  of  parts  I  and  II,  respectively,  cited  here  will 
be  preceded  by  the  numerals  I  and  II.  Given  a  plane  wave 
with  a  propagation  constant  k  *  2r/X  and  a  propagation  unit 
vector  fl„  Eq.  (1. 10a)  gives 

q(Z,  a,  at)  *  -({♦  -  fA)/fc  +  [sin  8,(Z,  cos  v,  +  Z,  sin 

-cos  «,]/(!  +  ZI2  +  ZV2)1--  (1) 


r 
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Fig.  1.  Geometry  of  the  problem. 
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<b)  (d) 

Fig.  2.  Results  of  Ref.  3,  Figs.  10-12,  reproduced,  using  our  computational  code,  with  a  perfectly  reflecting  paraboloidal  mirror  with  a  focal 
length  of  R! 2  •  50  cm  and  a  height  of  d  «  10  cm.  The  following  notation  is  used:  v  ■  kx  sin  o  on  the  jr  axis  and  u  “  ky  sin  a  on  they  axis,  sin  a  m 
(2Rd)'n/(R/2),  and  u  »  kz  sin- «.  (a)  Contours  of  the  reflected  energy  in  the  xy  plane  at  the  focus,  z  =  -R/2  (the  contour  labels  are  login  of  the 
reflected  energy  normalized  to  the  maximum  incident  energy),  (b)  Surface  plot  of  the  energy  in  the  xy  focal  plane,  (c),  (d)  Correspond  to  (a) 
and  (b),  respectively,  but  in  the  xz  plane  with  z  “  0  in  these  figures  corresponding  to  the  focal  point  in  Fig.  1. 


Equations  (1.8)  then  give  'HZ,  a,  w)  and  A(Z,  it,  w)  in  terms  of 
the  local  angle  of  incidence  defined  below  by  Eq.  (6). 

The  incident  electromagnetic  fields  have  field  strength 
distribution  E,(Z)  and  aberration  function  u(Z)  given  as 
follows. 

For  a  plane  wave, 

E,(Z)  =  F0  exp(-ikw)  (2a) 

and 

u  -  0.  (2b) 

For  a  Gaussian  beam,6 

E,(Z)  *  £0(ujn/u;,)exp{-(r/iz'1):!  -  ik(w  -  lul)]  (2c) 

and 


Here,  w  is  the  axial  distance  along  the  laser  beam,  measured 
from  the  beam  waist;  is  the  beam-waist  parameter;  w,  is  a 
beam  parameter  defined  by  w,2  «  m02(l  +  (2w/kwn2)2\\  r  is 
the  radial  distance  in  the  beam  frame  of  coordinates;  and  p  is 
the  wave-front  curvature  at  the  distance  w  and  is  given  by  p 
=  w  +  (ku>o2/2)2/w.  Thus  the  Fs  in  Eqs.  (1.14)  are  deter¬ 
mined  as  follows: 

For  a  plane  wave, 

F,  *  sin  8,  cos  (3a) 

F.,  *  sin  6,  sin  (3b) 

and 


u  »  [(l/fc)tan-l(2w/*eo02)  -  r1/2p\hr 


(2d) 


F,  =  cos  8,-. 

For  a  Gaussian  laser  beam, 


(3c) 
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F,  -  sin  0,  cos  *>,[1  +  (l/fc)tan"l(2u;Au;0->)  -  r/2(w  -  p)], 

(3d)  and 


Zs  =  -x/R 


(5a) 


F2  *  sin  9,  sin  *>,(1  +  (l/k)tan~l(2w/kw02)  -  r/2(w  -  p)], 

(3e) 

and 

F3  *  cos  0,(1  +  (l/fc)tan-1(2W)mv)  -  r72 (w  -  p)\.  (30 

Now,  the  coordinate  function  of  a  paraboloidal  mirror 
(Fig.  1)  can  be  written  in  the  beam  frame  (Z  »  —  z)  as 

Z  *  ~(xJ  +  y2)/2R  for  -d  <  Z  <  0,  (4) 

where  R/2  is  the  focal  length  of  the  paraboloidal  mirror. 
The  origin  of  the  frame  of  coordinates  is  taken  at  the  vertex 
of  the  paraboloid,  and  the  z  axis  is  its  axis  of  revolution.  The 
spatial  derivatives  of  Z  are 


Z,.  ■  —y/R.  (5b) 

The  local  angle  of  incidence  is  defined  by 

cos  0jn  =  (-sin  0,(Z,  cos  <e,  +  Zv  sin  <p,)  +  cos0,]/ 

(1  +  Z 2  +  Z2)m.  (6) 

Substituting  the  Fs,  Z%,  and  Z,  from  Eqs.  (3)  and  (5)  into 
Eqs.  (1)  and  (1.15)  gives  the  real  part  of  the  complex  refrac¬ 
tive  index,  ►,  for  tb»  surface  Z(x,y)  as  a  function  in  Z,  Z„  Z,, 
a,  and  u. 

The  Fresnel  coefficients  of  reflection  and  transmission, 
written  in  terms  of  v,  are  then  determined  in  a  form  similar 
to  Eqs.  (11.9).  Then  Eqs.  (1.19)  give  the  various  field  vectors 
e,‘,  h,\  ep't  and  hp‘,  and  Eqs.  (1.21)  and  (1.22)  give  the  local 
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F 


Fig.  3.  The  S's  [surface  currents  at  an  aluminum  paraboloidal 
mirror  surface,  Eqs.  (10)].  (a)-(c)  Are  for  0,  »  0  incidence, 

and  (d),  (e)  are  for  45*  incidence  in  the*r,  ”  0  plane.  (a),(d),S|.  (b),(e) 
S2.  (c)  Si. 


currents  on  the  surface  Z,  J,‘  and  Jp‘.  Thus  both  absorption 
and  reflecting  parameters  at  the  surface  are  fully  deter¬ 
mined. 

To  simplify  numerical  calculations,  only  the  case  of  a 
plane  wave  incident  at  0,  «  0,  *>,  *  0,  and  a  =  0  will  be 
considered.  For  this  case,  Eqs.  (1.18)  give 

t  *  [lZy  cos  9,  -  /(sin  0,  +  Zx  cos  9,)  -  £Z,  sin  0,]/C  sin  0in, 

(7a) 

p  -  |t]sin  0,(ZV2  +  1)  +  Zx  cos  0,]  +  /[— Z,Z,  sin  9,  +  Z}  cos  0(] 


+  £[cos  0,(ZI2  +  Z2)  +  Zx  sin  9JI/C2  sin  0in,  (7b) 

cos  0in  “  (cos  0,  —  Zx  sin  0,)/(l  +  Zx  +  Z2)m,  (7c) 

d  •  f  »  -ZJC  sin  0in,  (7d) 

A  •  A  *  (Zx  cos  0,  +  sin  0,)/C,  (7e) 

and 

d  •  P  *  [Z,( 2  sin2  0,  +  Zx  sin  0,  cos  0,  -  1) 

-  sin  0,  cos  0,]/C2  sin  0in.  (70 


Thus  all  the  field  components  in  the  incident  partial  electro¬ 
magnetic  waves  at  the  surface  Z(x,  y)  are  also  fully  deter¬ 
mined. 

At  a  point  P(xp,  yp,  zp)  external  to  the  surface  Z(x,  y),  the 
reflected  field  is  given  by  Eq.  (1.30).  The  unit  vectors  A  and 
Ap  are  given  by 

A  -  (~ZX!  -  Zj  +  *)/(l  +  Zx2  +  Z2)m  (8a) 

and 

Ap  =  [(x  -  x p)i  +  (y  -  yp)j  +  (Z  -  zp)&]/D,  (8b) 

where  D  -  ((x  -  xp)2  +  (y  -  y„)2  +  (Z  -  zp)2\,/2. 
Substituting  the  eJ’s  and  h'’s  from  Eqs.  (1.19)  into  Eqs. 


(1.31)  and  (1.32),  we  get  expressions  similar  to  Eqs.  (11.17)  for 
the  surface  currents  as 


A  X  e,.  =  (E,/CHSj  +  S.J  +  Sfi) 

(9a) 

and 

A  X  h,  =  {EJrjCHTj  +  Tj  +  Tfe). 

(9b) 

The  S’s  and  the  T s  (the  dimensionless  parameters  that 
account  for  the  physical  and  geometrical  properties  of  the 
surface,  the  degree  of  polarization,  and  the  radiation  wave¬ 
length)  are  now  defined  as 

S,  -  Wl-ZyLt  -  L2)  +  cos9,n(l  -  rtp)M„ 

(10a) 

S2  =  W(ZxL3  +  Lj)  +  cos0,„(l  —  Rp)M2, 

(10b) 

S3  -  W(-ZxL2  +  Z,X.)  +  cos 0jn(  1  -  Rp)M), 

(10c) 

T,  -  [(1  +  Rp)/C](-Z,M:X  -  M2)  -  ULX , 

(lOd) 

Ty  =  [(1  +  ffp)/Cl(Z,Af.t  +  A/,)  -  L'L.y, 

(10e> 

and 

T3  *  [(1  +  Rp)/C](-Z,M,  +  ZyMx)  -  UL). 

(lOf) 

By  executing  the  vectorial  multiplication  in  Eqs.  (1.31) 
and  (1.32),  we  get 

lV»d-t(l  +/?,)/C  sin  8m, 

(11a) 

t/  «  d  •  £(1  -  R,)cos  0in/sin  0jn. 

(lib) 

From  the  expression  of  e,',  Eq.  (I.19a),  we  get 

L,  *  Zv  cos  0,, 

(11c) 

L2  •  —(sin  0,  +  Zx  cos  0,), 

(lid) 

and 

L3  *  -Zv  sin  0,. 

(lie) 

From  the  expression  of  hp',  Eq.  (I.19d),  we  get 

M,  “  -K2  cos  0„ 

(110 

M2  =  K,  cos  0,  -  K3  sin  0,, 

dig) 

Fig.4.  Effect  of  th«  relative  central  obstruction  of  the  paraboloidal 
mirror  on  the  focal  intensity  (xp  “  0,  yp  »  0,  zp  “  R)  for  four  metals. 
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Fig.  5.  Intensity  distribution  along  the  x  axis  in  the  focal  plane  for 
the  four  metals. 


Fig.  6.  Intensity  distribution  along  the  z  axis  in  the  focal  volume 
for  the  four  metals,  z  «  0  in  this  figure  corresponds  to  the  focal 
point  in  Fig.  1. 


and 

A/3  =  K2  sin  0,.  (llh) 

From  the  expression  of  ep',  Eq.  (1. 19c),  we  get 
K ,  *  -&  •  hZx  +  &  •  /)[sin  9,{Z2  +  1)  +  Zx  cos  9J/C  sin  0in, 

(lii) 

K-2  =  -a  ■  hZ,  +  a  •  p(-ZxZy  sin  d,  +  Zy  cos  0,)/C  sin  0in, 

(llj) 

and 

K3  *  d  •  h  +  d  •  0[(Z,2  +  Z}2)cos  0,  +  Zx  sin  0J/C  sin  0in. 

(llk) 

Substituting  the  surface  derivatives  from  Eqs.  (5)  and  the 
vectorial  products  from  Eqs.  (7)  into  Eqs.  (10)  and  (11),  we 
get  expressions  for  the  S’s  and  the  T’s  in  terms  of  the  surface 
coordinates.  Then  the  surface  currents  that  contribute  to 


the  reflected  fields  at  P  are  given  by  Eqs.  (11.20).  The  three 
components  of  the  reflected  electric  field  at  P(xp,  yp,  zp)  are 
obtained  by  substituting  the  surface  currents  into  Eq.  (1.30). 
The  results  are 

E,X(P)  =  -[ifc  exp(-ikDu)/4irDu\  jj  lE,(Z)/D2jcos2  0,„ 

x  l[(y  -  y„)S:)  -  (Z  -  zp)S,]D  -  [(x  -  xp) 

X  (Z  -  zp)T,  -  [<y  -  yp)2  +  (Z  -  zp)2\Tx 
+  (x  -  xp)(y  -  yp)7'2]|exp(ife4>)dxdy,  (12a) 

E,y(P)  =  ~[ik  exp(—ikD0)/4*D0]  j  j  [E,(Z)/P2  |cos2  0ln 

X  |[(Z  -  zp)S,  -  (x  -  xp)S3]D  -  [(x  -  Xp) 
x  (y  -  yp)T,  -  [(x  -  xp)2  +  (Z  -  zp)2]T2 
+  (y  -  yp)(Z  -  xp)r3]|exp(i*«l>),  (12b) 

and 


Fig.  7.  Effect  of  the  relative  obstruction  on  the  z -axial  intensity 
distribution  for  an  aluminum  mirror,  z  m  0  in  this  figure  corre¬ 
sponds  to  the  focal  point  in  Fig.  1. 


Fig.  8.  Effect  of  the  relative  obstruction  on  the  t-axial  intensity 
distributions  for  an  aluminum  mirror. 
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E„(P)  =  -[/*  exp(-i7?D„)/4ir£>„)  j  j  {E,(Z)/D-]co s'2  0in 

X  |((x  -  xp)S.,  -  (y  -  >'P)S,]D  —  |Cv  —  >p) 

X  (Z  -  zp)T.,  -  [(x  -  Xp)~  +  (y-  yp)2]T:l 
+  (x  -  xp)(Z  -  xp)7’l]|exp(t^‘l>),  (12c) 

where  D()  is  the  distance  from  the  origin  to  the  observation 
point  and  4>  is  given  by 

+  =  r  •  rlp  =  (x(x  -  xp)  +  y(y  -  >p)  +  Z(Z  -  zp)]/D.  (12d) 

The  solid  angle  Q,  is  defined  in  Ref.  2. 

The  integrals  in  Eqs.  (12)  are  the  Stratton-Chu-Silver 
intergrals.  They  offer  the  following  advantages  over  the 
Kirchhoff  integrals:  (a)  the  physical  properties  of  the  re¬ 
flector  are  explicitly  included  in  the  reflecting  kernel,  while 
that  used  by  Barakat,  q(6)  =  2/(1  +  cos  0),  ignores  any 
surface  effect;  (b)  the  geometrical  surface  function  is  also 
explicitly  included  in  the  reflecting  kernel,  which  permits 
the  treatment  of  any  reflecting  geometry;  (c)  the  incident 
rays  need  not  be  parallel,  meridional,  or  symmetrical;  (d) 
integration  is  carried  out  on  the  surface  function,  not  on  a 
Gaussian  sphere,  permitting  the  treatment  of  asymmetrical 
focusing;  and  (e)  the  effect  of  the  thickness  of  the  metallic 
layer  can  easily  be  investigated  in  a  parametric  study  for  the 
engineering  design  of  paraboloidal  mirrors.  The  last-  named 
issue  was  treated  in  Ref.  7  by  using  the  Debye  potential 
solution  of  the  scalar  wave  equation.  The  Debye  potential 
approach  can  be  partly  applied  to  our  problem  if  the  thick¬ 
ness  of  the  metallic  coating  is  to  be  optimized.  In  the  region 
of  the  metallic  coating  of  few  tens  of  angstroms,  ray-tracing 
geometric  optics  is  of  no  use. 


the  focal  intensity  is  maximum.  Copper  and  silver,  as  ex¬ 
pected,  show  the  highest  focal  intensity,  while  aluminum 
shows  the  lowest  focal  intensity.  The  peaking  of  the  focal 
intensity  that  is  produced  by  using  a  specific  ring  aperture 
was  proved  by  others,  both  experimentally  and  theoretically, 
and  was  used  in  Ref.  8  to  produce  a  diffraction-free  beam. 

Figure  5  shows  the  symmetry  of  the  intensity  distribution 
along  the  x  and  z  axes  in  the  focal  volume  for  the  four  metals. 
Figure  6  shows  the  asymmetry  along  the  z-axis  intensity 
distribution  for  the  four  metals. 

Figures  7  and  8  show  the  effect  of  the  relative  obstruction 
on  the  axial  and  perpendicular  intensity  distributions  for  an 
aluminum  mirror. 

CONCLUSION 

In  an  application  of  the  theory  developed  in  part  I,  an  analy¬ 
sis  of  the  structure  of  the  electromagnetic  field  spectrum  at 
the  focal  volume  of  an  aluminum-coated  paraboloidal  mirror 
is  presented.  This  is  done  by  replacing  the  Kirchhoff  dif¬ 
fraction  integral  by  the  Stratton-Chu-Silver  integral.  The 
latter  is  more  appropriate  for  the  study  of  the  diffraction 
problems  by  reflectors  of  complex  shapes  and  physical  pa¬ 
rameters.  This  permits  the  inclusion  of  the  surface  physical 
and  geometrical  parameters  and  the  full  aberration  function 
in  the  reflecting  kernel.  The  numerical  results  obtained  are 
in  good  agreement  with  Barakat's  results  in  the  special  case 
of  a  perfectly  reflecting  mirror.  We  have  demonstrated  the 
effects  of  the  central  obstruction  on  the  focal  intensity  dis¬ 
tribution  for  four  practical  metallic  mirrors.  The  numerical 
results  demonstrate  the  use  of  the  analytical  expressions 
obtained  in  Eqs.  (12)  in  the  material  selection  and  design  of 
real  telescopic  mirrors  with  metallic  coating. 


NUMERICAL  RESULTS 

First,  the  results  obtained  by  Barakat  (Ref.  3,  Figs.  10-12) 
are  reproduced,  using  the  Kirchhoff  diffraction  method  with 
a  perfectly  reflecting  paraboloidal  mirror  with  a  focal  length 
of  R/2  *  50  cm  and  a  height  of  d  =  10  cm  along  the  mirror 
axis  of  revolution.*  The  Kirchhoff  integrals  were  evaluated 
by  a  numerical  double-integration  algorithm,  using  eight- 
point  Newton-Cotes  quadratures.  We  evaluated  the  inte¬ 
grals  without  using  the  Bessel  function  formulas  used  in 
Refs.  3  and  5.  Our  results  are  in  a  good  agreement  with 
those  obtained  in  Ref.  3  and  are  displayed  in  Fig.  2. 

In  Fig.  3  the  surface  currents,  the  S’s,  are  shown  for  both 
symmetrical  and  oblique  incidence.  The  figures  are  intend¬ 
ed  to  demonstrate  the  response  of  the  physical  parameters  of 
a  metallic  mirror  to  the  two  angles  of  incidence.  In  Fig.  3 
and  the  figures  that  follow,  A  *  1.315  jam. 

Figure  4  shows  the  effect  of  the  relative  central  obstruc¬ 
tion  of  the  paraboloidal  mirror  on  the  focal  intensity  (xp  =  0, 
yP  m0,zpm  R)  for  four  metals.  At  0.55,  central  obstruction 
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